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ABSTRACT 

Using  recently  published  results  of  intrinsic  and  free  carrier  nonlinear  absorption  coefficients  in  InP,  nonlinear 
refraction  was  investigated  at  1.064  ym  using  ns  duration  lasers  to  characterize  refraction  from  generated  free 
carriers.  A  phase  retrieval  algorithm  was  implemented  to  determine  the  amplitude  and  phase  profiles  of  the  inci¬ 
dent  beam.  Accurate  spatial  and  temporal  profiles  of  the  incident  held  were  used  to  model  nonlinear  propagation 
through  and  beyond  the  sample.  With  the  sample  held  fixed  at  focus  and  the  incident  energy  increased,  images 
of  the  transmitted  beam  a  fixed  distance  away  were  recorded  as  a  function  of  irradiance.  Excellent  agreement 
was  observed  between  recorded  beam  images  and  those  generated  from  the  numerical  model. 
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1.  INTRODUCTION 

Nonlinear  absorption  (NLA)  in  semiconductors  is  a  coupled  process  where  intrinsic  two  photon  absorption  (2PA) 
generates  free  carriers.  For  a  laser  of  sufficient  temporal  duration,  the  free  carriers  generated  locally  modify  both 
the  linear  absorption  and  refractive  index  within  the  material  while  the  laser  pulse  is  present.  Nonlinear  refraction 
(NLR)  due  to  free  carriers  has  long  been  reported1  and  values  for  InP  obtained  using  closed  aperture  Z-scans 
are  in  the  literature.2 

As  laser  pulse  durations  increase  from  sub  picosecond  to  nanoseconds  or  longer,  free  carrier  effects  become 
the  dominant  NLA  and  NLR  mechanisms.  Oftentimes  accurate  analysis  of  NLR  is  hampered  by  inaccurate 
knowledge  of  the  spatial  amplitude  and  phase  of  the  incident  beam.  This  has  been  overcome  by  use  of  a  phase 
retrieval  method.  By  capturing  images  of  the  incident  beam  as  well  as  a  second  image  located  a  known  distance 
away,  a  method  based  on  the  Gerchberg  Saxton3  algorithm  has  been  implemented  to  retrieve  the  spatial  phase 
profile  of  the  incident  beam.  With  the  amplitude,  phase  and  temporal  profile  of  the  incident  beam  known, 
the  held  was  numerically  modeled  through  the  sample  and  then  propagated  a  known  distance  beyond.  Both 
total  energy  transmission  and  spatial  beam  profiles  as  functions  of  incident  irradiance  were  recorded  and  good 
agreement  with  the  model  was  observed. 

2.  THEORY  AND  NUMERICAL  MODELING 
2.1  Coupled  Nonlinear  Beam  Propagation  and  Rate  Equations 

Equation  1  describes  the  propagation  of  the  complex  electro- magnetic  held  amplitude,  A(x,  y,  £),  through  a  third 
order  nonlinear  medium  along  the  £  direction  using  only  the  slowly  varying  amplitude  approximation. 

dA  i  _ n  A  2iir  (  A  dn  A  m  \  ,  A 

X  =  W0V±A  +  X  yre/A7V  +  dfAT  +  ^I)A~(a  +  a-»sAN  +  (31)-  (1) 

The  parameter  aref  is  similar  to  the  thermo-optic  coefficient,  and  relates  the  change  in  refractive  index 
due  to  generated  carriers,  The  instantaneous,  or  bound  electronic,  change  in  the  refractive  index  due  to 
irradiance,  /,  is  given  as  7.  While  included  for  completeness,  since  the  pulse  duration  studied  here  is  on  the 
order  of  nanoseconds,  NLR  due  to  7  is  much  weaker  than  free  carrier  refraction  effects  since  the  irradiances 
are  relatively  low.  One  and  two  photon  absorption  coefficients  are  a  and  (3  respectively  and  the  free  carrier 
absorption  cross  section  is  cra^s. 
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Generated  carrier  density,  AN(x,y,t),  and  temperature  change,  A T(x,y,t)  are  obtained  from  the  following 
rate  equations: 

dAN  _  olI  f3I2  AN 

dt  hv  +  2hv  r 

(^N~  =  —  [(a  +  aabsAN)I  +  f3I2] 

at  pc  L  J 

The  carrier  recombination  lifetime  r  is  much  longer  than  the  pulse  duration  used  here.  The  density  is  p  while  c 
denotes  the  specific  heat  of  the  medium.  The  2PA  coefficient,  /?,  and  the  free  carrier  absorption  cross  section, 
c Tabs  are  temperature  dependent  parameters.4  A  paper  detailing  the  theory  and  temperature  dependent  [3  and 
(Jabs  parameters  for  InP  is  in  preparation. 

2.2  Numerical  Method 

To  model  the  propagation  of  an  electro  magnetic  field  through  the  sample,  Eq.  (1)  is  numerically  solved  via 
the  operator  splitting  method.5  The  right  side  of  the  equation  is  separated  into  linear  and  nonlinear  operations 
which  are  solved  independently.  First  the  sample  length  is  divided  into  a  number  of  slices  of  length  A z.  Both  the 
amplitude  and  phase  of  the  incident  field  A(x,  y ,  t)  at  z ■.*==  Z\  are  known  from  the  experimental  conditions  and 
transmission  losses  are  accounted  for  to  determine  the  field  just  inside  the  sample.  Next,  with  the  nonlinearities 
turned  off,  diffraction  over  a  distance  4r  is  solved  for  each  time  slice  by  use  of  the  two  dimensional  fast  Fourier 
transform  beam  propagation  method  (FFT-BPM).6,7  Then  diffraction  is  set  to  zero  and  the  rate  equations  are 
solved  at  this  £  location  via  finite  difference  approximations.  Updated  values  for  A N(x,y,t)  and  A T(x,y,t)  are 
used  to  modify  the  electric  field.  Lastly  the  field  is  propagated  4r  and  the  process  is  repeated  until  the  pulse 
reaches  the  end  of  the  sample.  At  the  exit  of  the  sample  the  field  is  modified  by  the  transmission  coefficient  and 
is  temporally  and  spatially  integrated  to  calculate  the  output  energy.  Using  the  FFT-BPM  allows  propagation 
of  the  field  at  the  exit  of  the  sample  to  any  measurement  plane  after  the  sample.  Each  time  slice  of  the  field  is 
propagated  and  then  temporally  integrated  to  determine  a  fluence  distribution  for  comparison  with  the  camera 
image  recorded  at  the  same  plane. 

2.3  Determination  of  the  Incident  Amplitude  and  Phase 

To  determine  the  incident  spatial  field,  A(x,  y,  z  =  Zi,  £),  a  phase  retrieval  method  similar  to  Pandey  et  al.8  was 
used  but  with  the  addition  of  a  discrete  wavelet  transform  (DWT).  Beam  images  at  two  locations,  the  sample 
input  plane,  Zi,  and  the  detector  plane,  Z2,  were  recorded.  After  background  corrections,  the  square  root  of  the 
beam  images  are  taken  to  determine  the  field  amplitudes  at  each  plane.  An  initial  guess  is  taken  for  the  phase 
at  Z\  and  the  complex  amplitude  is  then  propagated  to  Z2.  The  phase  and  amplitude  are  separated  and  the 
resulting  propagated  amplitude  is  replaced  with  the  known  amplitude.  A  DWT  is  used  on  the  retrieved  phase 
to  filter  and  smooth  the  data  and  minimize  noise  effects.  A  complex  field  at  Z2  based  on  the  original  known 
amplitude  and  retrieved,  filtered  phase  is  back  propagated  to  Z\.  This  process  is  repeated  until  the  errors  are 
below  a  preset  threshold  value.  To  verify  the  accuracy  of  the  retrieved  phase,  the  field  is  propagated  to  a  third 
plane  beyond  Z2  and  compared  with  the  recorded  beam  image  there. 

The  incident  temporal  distribution  of  A(x,  y ,  t)  was  measured  using  a  fast  photodiode  and  oscilloscope 
(4  GHz).  The  pulse  shape  was  normalized  and  a  cubic  spline  used  to  decrease  the  number  of  sampling  points 
while  not  losing  any  temporal  features.  The  final  number  of  time  points  sets  the  number  of  time  slices  in  the 
simulated  amplitude. 

3.  RESULTS  AND  ANALYSIS 

3.1  Experimental  Layout 

The  experimental  set-up  shown  in  Fig.  1  consists  of  a  flashlamp  pumped,  electro-optically  q-switched  Nd:YAG 
laser  at  1.064  pm  (Quantel  Brilliant)  with  a  pulse  duration  of  3.7  ns  (FWHM).  No  spatial  filtering  of  the  beam 
was  done.  The  beam  was  sent  though  a  half  wave  plate  and  polarizer  for  attenuation.  The  half  wave  plate  was 
mounted  on  a  computer  controlled  rotation  stage.  Next  the  beam  was  reflected  by  a  45°  partial  reflector.  The 


(2) 

(3) 


2 


HWP&PoL  Cube 


InP  sample 
at  focus,  Z1 


Nd:YAG  Laser 


ChB 

Etrans  Monitor 

tfi - 1 


2X  telescope 


□ 


Camera  atZ 


CaF2  wedges  (x2) 


Figure  1.  Experimental  layout 


transmitted  beam  was  sent  to  a  pyroelectric  detector  used  to  monitor  the  incident  energy  (Ch.  A  detector  in 
Fig.  1).  Next  a  2X  down  telescope  was  used  to  increase  the  irradiance  at  the  sample.  The  beam  was  loosely 
collimated  after  the  down  telescope  to  ensure  that  the  sample  thickness  was  thin  compared  to  the  diffraction 
distance  of  the  beam.  The  irradiance  scan  method9  was  used  where  the  sample  was  kept  fixed  at  focus  while  the 
incident  energy  was  varied. 

A  pair  of  CaF2  wedges  were  placed  immediately  after  the  sample.  The  front  surface  reflections  from  the 
wedges  were  sent  to  a  camera  located  at  Z2  while  the  transmitted  beam  energy  through  the  first  wedge  was 
measured  by  a  second  pyroelectric  energy  detector  (Ch.  B).  The  distance  between  planes  Zi  and  Z2  was  650  mm 
with  Z2  located  within  the  Rayleigh  range  of  the  beam.  Placing  the  detector  in  the  near  field  maximizes  the 
amount  of  observed  signal  due  to  NLR.  If  the  detector  is  placed  in  the  far  field,  less  NLR  effects  are  observed 
since  linear  diffraction  due  to  free  space  propagation  will  dominate. 

3.2  Nonlinear  Absorption  Results 

Figure  2  shows  the  nonlinear  transmission  measured  as  the  ratio  of  detectors  B/A  as  a  function  of  incident 
energy. 


Einc(J) 

Figure  2.  Nonlinear  absorption  in  InP  at  1.064  /xm  using  nsec  pulses. 


In  Fig.  2  the  circles  are  experimental  data,  dashed  line  is  numerical  modeling  results  using  only  2PA  and  solid 
line  is  modeling  results  including  both  2PA  +  FCA.  The  InP  sample  was  1.0  mm  thick  and  was  anti-reflection 
coated.  The  incident  beam  was  elliptical  and  approximately  Gaussian.  Fitting  a  Gaussian  to  the  beam  irradiance 
gives  1/e  radii  of  1.0  and  0.5  mm  in  the  horizontal  and  vertical  directions.  The  pulse  width  was  2.8  nsec  FWHM. 
The  values  of  (3  =  25.5  cm/GW  and  aabs  =  1.5  x  10-17  cm2  were  used  in  the  numerical  modeling.10 
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3.3  Nonlinear  Refraction  Results 

Measurement  of  nonlinear  refraction  was  done  by  recording  camera  images  at  plane  Z 2  as  a  function  of  the 
incident  irradiance.  Figure  3  shows  the  beam  recorded  at  Z2  (left)  and  the  result  of  the  numerical  model  (center) 


Figure  3.  Beam  images  from  camera  (left)  and  modeling  results  (center)  at  plane  Z2  for  incident  irradiance  of  79  MW/cm2. 
Right  column  shows  a  comparison  of  the  horizontal  line  profiles  through  the  centers  of  both  images. 

for  an  incident  irradiance  of  79  MW/cm2.  The  graph  (right)  shows  a  horizontal  line  profile  through  the  center 
of  both  images.  The  free  carrier  refraction  coefficient  was  determined  by  iteratively  adjusting  the  value  of  the 
parameter  while  minimizing  the  errors  between  the  recorded  images  and  numerical  result  over  all  irradiances. 
Note  that  the  beam  images  from  the  camera  and  model  are  not  on  the  same  intensity  scale. 


4.  CONCLUSION 

With  recently  published  results  of  intrinsic  and  free  carrier  nonlinear  absorption  coefficients  in  InP,  nonlinear 
refraction  was  investigated  at  1.064  /am  using  ns  duration  lasers  in  order  to  characterize  refraction  from  generated 
free  carriers.  A  phase  retrieval  algorithm  was  implemented  to  determine  the  amplitude  and  phase  profiles  of 
the  incident  beams.  Accurate  spatial  and  temporal  profiles  of  the  incident  beams  were  used  to  model  nonlinear 
propagation  through  the  sample  and  linear  propagation  beyond.  With  the  sample  held  fixed  at  focus,  images 
of  the  transmitted  beam  a  fixed  distance  away  were  recorded  as  a  function  of  incident  irradiance.  Excellent 
agreement  was  observed  between  recorded  beam  images  and  those  generated  from  the  numerical  model. 
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